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We adopt the Dirac model for graphene and calculate the Casimir interaction 
energy between a plane suspended graphene sample and a parallel plane ideal 
conductor. We employ both the Quantum Field Theory (QFT) approach, and 
the Lifshitz formula generalizations. The first approach turns out to be the 
leading order in the coupling constant of the second one. The Casimir interac- 
tion for this system appears to be rather weak but experimentally measurable. 
It exhibits a strong dependence on the mass of the quasi-particles in graphene. 
Present article is based on joint works. 
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1. Introduction 

Graphene is a (quasi) two dimensional hexagonal lattice of carbon atoms. 
It belongs to the most interesting materials in solid state physics now due 
to its exceptional properties and importance for nano technology.'^''* Here 
we consider the Casimir interaction between suspended graphene plane and 
parallel ideal conductor. This setup was considered in^"'' using a hydrody- 
namical model for the electrons in graphene following. ^'^ Later it became 
clear that this model does not describe the electronic properties specific to 
this novel material. 

Here we use a realistic and well-tested model where the quasi-particles in 
graphene are considered to be two-component Dirac fermions. This model 
incorporates the most essential and well-established properties of the their 
dynamics: the symmetries of the hexagonal lattice, the linearity of the spec- 
trum, a very small mass gap (if any), and a characteristic propagation ve- 
locity which is 1/300 of the speed of light. By construction, this model 
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should work below the energy scale of about leV, but even above this limit 
the optical properties of graphene are reproduced with a high precision. 
The action of the model, therefore, is given by 

St,= j (f x4>{l\idi - eAi) - m)i/., / = 0, 1, 2 (1) 

where 7' are just rescaled 2^2 gamma matrices, 7° = 7", 7^'^ = Ui?7^'^, 
7q = —(7*)^ = 1- In our units, h = and Fermi velocity vp — (300)^^. 

The value of the mass gap parameter m and mechanisms of its generation 
are under discussion. ^^"^"^ The upper limit on m is about Q.leV at most. 

The propagation of photons in the ambient 3 + 1 dimensional space is 
described by the Maxwell action 

Sm^ -^J d'^xFf^.F'''', Ai,z/-0,1,2,3. (2) 

In the following we shall suppose that the graphene sample occupies the 
plane = a > 0, and the conductor corresponds to a;^ = 0. 



2. QFT approach 

In the framework of QFT one evaluates the effective action in a theory 
described by the classical action 5d + Sm- Then the Casimir energy density 
per unit area of the surfaces at the leading order in the fine structure 
constant a is given by 

= _J_^^, (3) 

where T is time interval, and S is the area of the surface. The solid line 
denotes the fermion propagator in 2 + 1 dimensions (i.e., inside the graphene 
sample) , and the wavy line is the photon propagator in the ambient 3+1 
dimensional space subject to the perfect conductor boundary conditions at 

X^ = 0: ^0^3=0 = ^lU3=0 = ^21x3=0 = 93^31^3=0 = 0. 

The fermion loop in 2 + 1 dimensions has already been calculated in a 
number of papers. ^^^^^'"'^^ It gives the quadratic order in A of the effective 
action for electromagnetic field 

SMA) = A w,^^^^^ A = ^J ^A,{p)W'ip)Aip), (4) 

where 

n- = =^<(."-f ).r, (5) 
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is the polarization tensor in the lowest, one loop, order. Here 77™ = 
dia.g{l, vf,vf), p denotes the rescaled momenta pj = rjjPk, The func- 
tion i>(p) is model dependent, and for graphene it reads = 
2 (2mp — {p^ + 4m^)arctanh(p/2TO)) /p. We assume here that all possible 
parity-odd parts are canceled in 11. Possible effects invoked by their presence 
are considered in.^ 

To calculate the diagram (3) we only need to couple the kernel (5) 
to the photon propagator subject to conducting boundary conditions. In 
Fourier representation and for the Euclidean 3- momenta, i.e., after the Wick 
rotation p ^ pe = {ipa,Pij P2), the a-dependent part of the energy reads 

_1 f d'pE mPE) f d'pE "(Pf -2ap„ 

4 J (27r)3 J (27r)3 4p||pf 

(6) 

where we expanded Hj{pe) explicitly with help of (5), and = \pe\- 
3. Lifshitz formula approach 

One can also consider the system as described by effective theory of the 
electromagnetic field with the action + Sgg subject to the conduct- 
ing boundary conditions at a;^ = 0. Then at the surface of graphene, the 
Maxwell equations receive a singular contribution 

^|,F^"' + 5{x^ - a)Il''PAp = (7) 

following from SeS- Here we set H^^ = 11''^ = 0. This contribution is 
equivalent to imposing the matching conditions 

assuming that is continuous at x'^ = a. Now, one can forget the origin 
of and quantize (at least formally) the electromagnetic field subject to 
the conditions (8) at cc^ = a and to the conducting conditions at x"^ = 0. 

The original Lifshitz approach-'^^ was generalized^^'^^ for the interactions 
between two plane parallel interfaces separated by the distance a and pos- 
sessing arbitrary reflection coefficients rj,^ , r!^^ of the TE and TM 
electromagnetic modes on each of the surfaces 

= / - - ^-'''''^'Tlr^Tl)]- (9) 

For graphene with help of matching conditions (8) we can obtain at the 
Euclidean momenta 

^^^"2^,11+^$' '^™- 2pf + ap||$' ^'^> 
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while for the perfect conductor one has = —1, r^2i = 1- It is clear, that 
$ must be rotated to Euclidean momenta as well. Wc also note that the 
perfect conductor case is recovered from (10) in the formal limit $ oo. 

One can show by a direct computation that the energy £i, Eq. (6), coin- 
cides with the leading order in a perturbativc expansion of the Lifshitz 
formula (9)-(10), so that the two approaches are consistent. 

4. Results and discussion 

The formulae (6) and (9)-(10) are suitable for the numerical and asymptot- 
ical evaluation. The asymptotic expansion for short and long distances are 
readily obtained through uniform expansion of the integrand of (6), (9) 

Note that the asymptotics at large separations is of the first order in a 
while for small separations, it contains all powers of a through h{a,VF), 
for the real values of parameters in graphene /i(l/137, 1/300) w 0.024. 

Therefore wc sec that at large separations Casimir energy docs not turn 
into the ideal conductor case, while at small separation this case is indeed 
recovered. This is counter-intuitive since the main contribution at short 
separations shall come from the high momenta for which one would expect 
the graphene film to become transparent. We must also stress that this 
behavior is drastically different from that in the hydrodynamic model. ^-'^ 

For numerical evaluation we normalize the results to the Casimir energy 
£c = — a'-^ ^'^^ plane ideal conductors separated by the same distance 
a. The results of calculations are depicted at Fig. 1. The scale is defined by 
the mass parameter m. For m of the order of next nearest-neighbor hopping 
energy t', i.e., m = O.le^,^ ma = 1 corresponds to a = 1.97 micrometer. 

Thus, we can sec that the magnitude of the considered Casimir inter- 
action of graphene with a perfect conductor is rather small. Actual mea- 
surement of such weak forces is a challenging, but by no means hopeless, 
experimental problem. Strong dependence on the mass parameter m at large 
separation is also a characteristic feature of the Casimir force. Getting an 
independent measurement of m may be very important for our understand- 
ing of the electronic properties of graphene. The mass of quasi-particles in 
graphene is, probably, very tiny, which improves the detectability of the 
Casimir interaction since the energy increases with decreasing m. 
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Fig. 1. The relative Casimir energy densities Si/£c (solid line) and £i^/Sc (dashed 
line) as functions of ma. Insert shows a zoom of the small-distances region. 
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